Scaling Properties of Weak Chaos in Nonlinear Disordered Lattices 
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The Discrete Nonlinear Schroedinger Equation with a random potential in one dimension is stud- 
ied as a dynamical system. It is characterized by the length, the strength of the random potential 
and by the field density that determines the eflect of nonlinearity. The probability of the system to 
be regular is established numerically and found to be a scaling function. This property is used to 
calculate the asymptotic properties of the system in regimes beyond our computational power. 
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The Nonlinear Schrodinger Equation (NLSE), in ab- 
sence of a potential, in the continuum is an integrable 
problem |l[. It is relevant for the description of Bose- 
Einstein condensates (where it is known as the Gross- 
Pitaevskii Equation [2]) as well as for classical nonlin- 
ear optics and plasmas [3!|. For the linear Schrodinger 
Equation with a random potential in one dimension all 
the states are localized [J], as a manifestation of An- 
derson localization [5|. It is natural to ask what is the 
asymptotic behavior for the NLSE with a random po- 
tential, namely what is the outcome of the competition 
between Anderson localization and nonlinearity. This 
question is directly related to several experimental sit- 
uations. For Bose-Einstein Condensates the nonlinear 
term in the NLSE models the interaction between the 
atoms while the random potential may be generated by 
lasers as in recent experiments \Q\. In nonlinear optics 
the NLSE adequately describes light propagation in a 
nonlinear medium where the randomness can be written 
with light that passes trough a diffuser 7] or by manu- 
facturing [8|. 

Here the exploration will be performed in the frame- 
work of a discrete one dimensional model (defined by 
([1])). The specific elementary question in this field is: 
will an initially localized wave packet spread to infin- 
ity in the long time limit? Extensive numerical simu- 
lations [Sf (where although the full control of errors is 
impossible, one argues that the average results are statis- 
tically meaningful) exhibit sub-diffusion, with the width 
of the wave packet growing in time as t" with a ~ 1/6. 
On the other hand, recently it was argued that eventually 
the spreading should stop and the dynamics is eventually 
almost- per iodic on a kind of KAM torus [lOJ. Rigorous 
studies 
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lead to the conjecture that in the strong dis- 
order limit the spreading is at most logarithmic in time, 
excluding sub-diffusion as the asymptotic behavior. Non- 
rigorous results based on perturbation theory extend this 
conjecture beyond the regime of strong disorder with the 
help of a bound on the remainder term of the pertur- 
bation series [12|, but the times available here at orders 



calculated so far turned out to be short compared with 
numerical calculations where sub-diffusion was found. A 
major difficulty in the exploration of this problem is that 
we do not know how far in space and time one should 
go so that the result can be considered asymptotic. One 
reason this problem is complicated is the fact that dur- 
ing the spreading the effective number of degrees of free- 
dom increases enhancing chaos, but their amplitude de- 
creases suppressing chaos. This motivated the present 
work which is designed to decide which of these com- 
peting effects wins. To address this issue we develop 
here a scaling theory of weak chaos in disordered non- 
linear lattices, expecting it will be useful for extending 
results beyond our computational ability. Scaling ap- 
proaches proved to be extremely powerful in equilibrium 
and nonequilibrium statistical physics ^] , and have been 
also very successful in understanding Anderson localiza- 
tion '1411. In this letter a scaling theory for the probability 
distribution to observe chaos or regularity, based on the 
computation of the largest Lyapunov exponent is devel- 
oped and tested for relatively small systems which are 
within our numerical power. Such a theory is expected 
to have predictive power when extended to infinite size. 
We study a nonlinear disordered medium described by 
the Discrete Anderson Nonlinear Schrodinger Equation 
(DANSE) model for a complex field 4'n{i)'- 



dt 



iTpn + Jii^n+1 + i'n-l) + \ Tpn \ V'ti 



(1) 



The hopping is J = (1 -f W)~^ while e„ are independent 
identically uniformly distributed in {—JW,JW). With 
this rescaling the eigenvalues E of the linear part satisfy 
l^-l < 1 + (f + W)^^. Hence for strong disorder (large W, 
that is the focus of our present work), the energies of the 
corresponding linear equation are practically in the inter- 
val (—1,1). Measuring the length scale of the eigenstates 
of the linear problem /i by the inverse participation num- 
ber /i^^ = J2k 1*^*1*' "^"2 find that /x « H- M^~^ for large 
W. By scaling the amplitude of the field we set the co- 
efficient of the nonlinear term in ([1]) to one. We consider 



DANSE model on a lattice (ring) of length L, with peri- 
odic boundary conditions. While DANSE model has two 
integrals of motion, the norm of the field N = ^^ |V'n|^ 
and the total energy, only the norm is important here; 
and in our treatment we do not control the energy which 
is always chosen to be close to zero (center of the band 
of the corresponding linear system). 

Our goal in this letter is to study properties of the 
dynamics as a function of all relevant parameters: disor- 
der W, norm N, and lattice length L. Equivalently, we 
introduce the density p = N/L and consider the depen- 
dence on two intensive parameters W and p and on the 
extensive parameter L. 

We characterize the dynamics by means of the largest 
Lyapunov exponent. For a particular realization of dis- 
order {e„}, we followed a dynamical solution of Eq. ([T]) 
starting with a uniform initial field, and for this trajec- 
tory calculated the largest Lyapunov exponent A by a 
standard method. This calculation was repeated for a 
large number of realizations of disorder, and the parame- 
ters above. As a result, we can construct a distribution of 
largest Lyapunov exponents over the realizations of dis- 
order, for given macroscopic parameters W,p,L. Several 
examples of these distributions are presented in Fig. [1] 
We can see here that for very small densities p, where the 
nonlinearity in DANSE is very small, all the Lyapunov 
exponents are small and close to ~ 10^^. For a regular 
dynamics the largest Lyapunov exponent should be ex- 
actly zero, but numerically, with a finite integrating time 
{Tmax = 10^ in our numerical simulation; control runs 
with Tjnax — 10^ showed no significant difference) such 
a small value essentially indicates regular (quasiperiodic) 
dynamics of the field. Contrary to this, for large densities 
we observe Lyapunov exponents in the range 0.01 — 0.1, 
what indicates for chaotic dynamics. For intermediate 
densities we see that for some realizations of disorder the 
dynamics is regular (A « 10~^), while for other realiza- 
tions much larger exponents are observed indicating weak 
chaos. As we want to perform a statistical analysis rather 
than going into details of particular dynamics for partic- 
ular realizations of disorder, we adopt the following op- 
erational definition to distinguish between regularity and 
chaos: we attribute all runs with A<5-10~^as regular 
ones, and all runs with A > 5 • 10~® as chaotic ones. In 
this way we define the quantity of our main interest in 
this letter, a probability to observe regular dynamics P. 

As explained above, we have determined the probabil- 
ity of regular dynamics P(p, W, L) as a function of pa- 
rameters of the model. Typical profiles of P for fixed 
disorder W and different values of p and L are depicted 
in Fig. [21 One can see a typical sigmoidal function with 
limits P — > 1 for p ^- and P — ;► for p — > cxi. 

We first concentrate on the dependence on the exten- 
sive parameter L. As one can see from Fig. [21 for the 
fixed W and p, the probability to observe regularity de- 
creases with the length L. This behavior can be under- 
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FIG. 1. Distributions of the largest Lyapunov exponent for 
L = 16 and W — 10, and three values of the field density. For 
p = 0.0003 all the realizations are regular, for p — 0.1 all are 
chaotic, but for p = 0.01 a part of realizations are chaotic. 
Vertical dashed line indicates the border A = 5 • 10~®. 
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FIG. 2. The values of P{p, W, L) for W = 10 and different L, 
as functions of p. For the same density p the probability to 
observe regularity decreases with L. 



stood as follows. Suppose we fix W, p and consider a 
lattice of large length L. Let us divide this lattice into 
(still large) subsystems of lengths Lq. How the probabil- 
ity to observe regularity on the large lattice P(p, W, L) 
is related to the corresponding probabilities for smaller 
lattices P(p, VF,Po)? It is reasonable to assume that to 
observe regularity in the whole lattice we need to have 
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FIG. 3. Rescaled according to Q probabilities to observe 
regular behavior, for different lattice sizes, (a): fixed disorder 
W = 10, dependence on density (the same data as in Fig. [2]); 
(b): fixed density p = 0.01, dependence on the disorder. 



all the subsystems regular, because any one chaotic sub- 
system will destroy regularity. This immediately leads us 
to the relation 



P{p,W,L)^[P{p,W,Lo)] 



L/Lo 



(2) 



Equation ^ implicitly assumes that chaos appears not 
due to an interaction between the subsystems, but in 
each subsystem (of length Lq) separately. This appears 
reasonable if the interaction between the subsystems is 
small, i.e. if their lengths are large compared to the 
length scale associated with localization in the linear 
problem: Lq ^ M- On these scales the various subsys- 
tems are statistically independent. This is the content 
of ([2]). It motivates the definition of the L-independent 
quantity: 



R{p,W)^[P{p,W,L)] 



l/L 



(3) 



We check the scaling relation (I2I3P in Figs. [3] and see that 
the data for lattices of sizes 16 < L < 128 collapse, so 
that R is independent of L. Remarkably, a short lattice 
with L = 8 obeys the scaling for large disorders but devi- 
ates significantly for small disorders M^ ^ 10; this corre- 
sponds to the expected validity condition that L should 
be larger than p, (the spatial size of eigenfunctions). 

The scaling relation ([2]) describes dependence on the 
extensive parameter L (and will allow us to extrapo- 
late results to long lattices beyond our numerical re- 
sources), so we can concentrate on considering dependen- 
cies on intensive parameters W and p. Therefore below 
we fix L = Lq = 16 and study the scaling properties of 
Po{p,W) — P{p,W,Lo). As this quantity is roughly a 
sigmoidal function of p (see Fig. |4^), it is convenient to 
perform a transformation to a new quantity Q{p, W) as 
Q = ^^°p . In this representation 



Pn = 



1 



Q{p,W) 

l + Qip,W) l + Q-^p,W) 



(4) 



so that the asymptotic behaviors Pq ^ 1 as Q ^ oo and 
Po — ?> as Q — ?> can be easily visualized as in Fig. |4)d. 
The next crucial observation is that the function 
Q{p,W) is not an arbitrary function of p and W, but 
it can be written in a scaling form 



Q 



w^^yw^j 



(5) 



where q(x) is as usual a singular function at its limits 
q{x) ^ cix~'' for small x, while q{x) ~ C2X~'' for large 
X. These are found from the straight lines in Fig. \5\ 
The data of Fig. |4] collapses to one curve as shown in 
Fig. [5] This is the numerical justification for ([5]). It also 
provides the values of the exponents a — (3 — 1.75, ( » 
9/4 = 2.25, T] « 5.2, a « 2.5 • lO"'^ and ca « 1.8 • IQ-^^ 
The existence of the scaling function enables us to an- 
alyze the behavior at various limits. Most interesting is 
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FIG. 4. Profiles of Po{p, W) vs p for different values of W (a) 
and the same data in terms of Q{p, W) (b). 




p/W 



FIG. 5. (Color online) The same data as in Fig. |4] but in 
scaled coordinates. The black dashed and red dotted lines, 
showing asymptotics for small and large arguments of q, have 
slopes (, = 9/4 and 77 = 5.2. 



the limit of small densities where Q is large (see Fig. [S]) 
and probability of regular behavior is close to one; corre- 
spondingly, the probability Pch to observe chaos is small. 
Using the general relation following from (|2I4I5P 



P{p,W,L) = 



W"q-' f-J^ 



we obtain 



P 



ch 



lOgP: 



LL^'W 






Using the asymptotics q{x) ^ cix ^ we get[l5| 
P,h ~ LLo^W^^'-^^p^c^^ 



(6) 



Now let us assume that we consider the states with 
the same fixed norm N on lattices of different length L. 
Then p = N/L and from Eq. ([6|) it follows 



P 



L^-^N^W^^^-^^ L-5/4^9/4 



ch 



Clio 



CiLoVF35/16 



(7) 



This quantity, as expected, grows with the norm N and 
decreases with the disorder W. We see that because 
C > 1, probabihty to observe chaos in large lattices at 
fixed norm tends to zero. This result may have impli- 
cations for the problem of spreading of an initially local 
wave packet in large lattices. In this setup the norm 
of the field is conserved, and the effective density de- 
creases in course of the spreading. If one assumes that 
the dynamics follows the scaling above (although we es- 
tablished it for a special setup of lattices of finite lengths 
with periodic boundary conditions), and if one assumes 
that chaos is essential for spreading, then one concludes 
that the spreading should eventually stop as the proba- 
bility to observe chaos eventually vanishes. To estimate, 
according to arguments above, at which length we can 
expect the spreading to stop, we have to start with large 
densities (where the probability to observe regular dy- 
namics is negligible) and to estimate, at what lattice size 
Lmax chaos extincts. Assume that this happens when 
the probability to observe chaos reaches some small level 
D. Substituting Pch ~ D in Eq. ([7]), we obtain for the 
following estimate for Lmax'- 



Assuming DLq 
we obtain 



N<^W-"{DLoCi)-- 



1 and substituting the constants found. 



,{N, W)k2-10^ ■ Ar9/5M^"'^/4 



(8) 



This is our estimation for the maximal spreading of a 
wave packet of norm A^ in a lattice with disorder W. 
Taking values typical for the numerical experiments [9|, 
namely A^ = 1, 1 <W < 10, we obtain from ([5]) L^ax in 
the range from 3 ■ 10* to 2 ■ 10^. This explains, why in the 
numerical experiments, where typically values L « 100 
are achieved, no saturation of the power law spreading is 
observed. 

In conclusion, we established a full scaling theory for 
the probability to observe chaos in disordered nonlinear 
Schrodinger lattice: the scaling with the extensive pa- 
rameter - lattice length - is given by ([3]), while the scaling 
dependence on intensive parameters (disorder and den- 
sity) is given by Eq. ([5]). The found scaling indices in- 
dicate, that for long lattices with the same norm chaos 
extincts and regularity prevails. Furthermore, we use the 
system presented here as a model for the chaotic region 
of high density typically observed at the initial stage of 
evolution in all numerical experiments [9|. In this con- 
text the scaling with the length allows us to estimate the 
maximal length that could be reached by spreading from 
initially local chaotic wave packet. Remarkably, the scal- 
ing relations established include large/small constants, 
what explains previous observations of energy spreading 
over extremely large time scales. Nevertheless, by apply- 



ing the scaling we were able to estimate the final stages 
of evolution from the studies of relatively small lattices. 
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